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In this article we combine the ideas introduced by us earher in various proportions to ar- 
rive at a simple and yet powerful means of studying single-particle properties of homogeneous 
Fermi systems in detail without making assumptions regarding the validity or otherwise of 
Fermi-liquid theory. Novelties include the exact forms of the momentum distribution and 
spectral functions in the intermediate density regime and a dielectric function that is sensi- 
tive to significant qualitative changes in single-particle properties. This time, care is taken 
to include the particle-hole mode by interpreting the sum over modes as a weighted sum 
with the dynamical structure factor appearing as the weight. 



I. INTRODUCTION 



■ Both the Hartree-FockH and the Random Phase approximations (RPA)a are widely used in the physics 
hterature to study Fermi systems. The Bogoliubov theorylj is the analog of the randontphase approxi- 

■ mation for Bose systems. This latter fact has been demonstrated in detail in our articleQ. The analysis 
presented here could be repeated for Bose systems as well. But let us focus on the more important Fermi 

. system. Here we try and address the question of validity of the Hartree-Fock approximation and the RPA. 

' This is important since we showed in our earlier work Ihat the simplest and most natural application of 

. bosonization was in the regime where the RPA was exacta. It is therefore quite pertinent to ask whether this 

approximation is adequate. It is well-known that there are many problems with the RPA. Iii particular, the 
J2 \ pair-correlation function has unphysical behaviour at short distances (it becomes negative )u. Therefore it is 

necessary to find better approximations. Many attempts have been made in the literature in this regard. The 
Hubbard approximation', the Singwi-Sjolander approach^ and many others discussed in the text by MahanEl 
are among the more prominent. Many of them have found to resolve this problem of short distance behavior 
Q ' of the pair correlation function. However it seems that none of them are able to address the question of 

. single-particle properties. This issue has acquired an urgency in the recent past given that a segment of the 

high Tc community is convinced that the non-superconducting state of these materials is non-Fermi liquid. 
Therefore being able to find a theory that describes the non-Fermi liquid regime accurately is desirable. In 
' this attempt, we shall adopt a two-track approach. On the one hand we reinterpret the RPA so that it be- 

Jh \ comes more widely applicable. We argue that while the mean-field approximation carried out on the density 

operator is the Hartree-Fock approximation, the RPA manifests itself as mean-field theory applied not to 
the density operator but to the number operator. The density operator measures how the electrons are dis- 
tributed in real space whereas the number operator measures how the electrons are distributed in momentum 
space. Just as the Hartree-Fock approximation is valid when fluctuations in the density of electrons at each 
point in real space is small, the RPA is valid when fluctuations in the momentum distribution of electrons is 
small compared with the average momentum distribution which measures the probability of an electron to 
possess a given momentum. We also introduce another mathematical transformation within the RPA-scheme 
and find that this maps a purely repulsive interaction to a purely attractive one but one that makes the 
particle exchange a momentum different from the usual one involved in two-body collisions. This unusual 
transformation allows us to define and compute yet another dielectric function of this system. A physical 
interpretation of this curious change in sign is given. A suitable combination of the purely repulsivc(both 
actually and apparently) and purely attractive(only apparently, the electrons still repel one another) is shown 
to be better than either one spearately. In the end, we point out how a motivated researcher may be able 
to compute the phase diagram of the homogeneous electron gas by appropriately combining these three 
ingredients, although we ourselves restrict our attention to simple analytically solvable cases. 

(i) Sea-bosons, indispensable for computing single-particle properties. 
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(ii) Generalised RPA and beyond (with fluctuations in the momentum distribution). Many new dielec- 
tric functions may be found here. 

(iii) Repulsion-attraction duality : the enigmatic transformation that allows us to view the purely repulsive 
interaction between electrons as being attractive and repulsive at the same time when viewed in the sea- 
boson language. 

In what follows, we try and combine the bosonization approach of our earlier articles with this notion 
of generalised RPA and repulsion-attraction duality and try and compute the single-particle properties. Just 
in order to whet the reader's appetite for the calculation of single-particle properties using sea-bosons, we 
pause here to explore another option, one that suggests itself quite naturally when one views RPA as a 
mean-field idea applied to the number operator rather than the density operator. 



II. SINGLE-PARTICLE PROPERTIES WITHOUT USING BOSONIZATION 

Let us try and write down some representativCpexamples when RPA and Hartree-Fock approximations are 
used. Let us take for example, the jellium modelQ. 

k q#0 

If we apply the mean-field idea on the density, that is, replace pq by (pq) then we get a hamiltonian that 
does not involve any coulomb interaction at all(apart from an additive constant). 

k q#0 



k q#0 

Therefore this approximation is bad in the extreme. However the random-phase approximation is still valid 
for this system. It is the mean-field idea applied not to the density but to the number operator. For this we 
have to rewrite the full hamiltonian given below, 

^k4ck + E E 4+q/2'^l'-q/2V+q/2Ck-q/2 (3) 

k q^^O k,k' 



This has to be replaced by, 



^ = ^0 - E 1^ E 4+q/2Ck'+q/2cI'_q/2Ck-q/2 (4) 



where 



Ho = Yl ^kclck - E ^ E "k+q/2"k-q/2 (5) 
k q#0 k 

where ~ Cj{j.Ck. Let us write, 

nk = ("-k) + Snu (6) 

We plug the above decomposition into Hq and find that if we neglect terms quadratic in the fiuctuations, 
we get what we are after, namely the generalised RPA called for simplicity as just RPA. 

Hrpa -E'o + E ^k^Ck (7) 
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The average occupation is 



ek^ e]^-^y{nk^q) (8) 



The chemical potential has to be fixed by making sure that, 

E("k) = {N) (10) 

k 

At zero temperature ^ ^ ep, this quantity is equal to the usual Fermi energy when Vq = 0. 



(rik) = eiep - e~k) = Oiep - fk + J] y K-q» (11) 

and is the Heaviside step function. |-| We can now demonstrate that the generalised RPA dielectric func- 
tion(the Lindhard dielectric functioiiD being a weak coupling limit of this) may be recovered using the 
following procedure. If one considers an extremely weak external perturbation applied to the system and 
follows the discussion in MahanO one arrives at the following formula for the dielectric function, 

^ ^ 1 , '^q ('^k+q/2> - ("k-q/2) 

^ ^ ^ <^k+q/2 + ek-q/2 

The only point to bear in mind is that we have to use the full interacting momentum distribution rather than 
just its noninteracting value. When one includes fluctuations in the momentum distribution however, the 
answer given later, is very different from the usual RPA or even the generalised R-PAi- This feature of having 
the full momentum distribution in the numerator may be found in our earlier worlafl. Thus we can see that 
there is a whole new set of approximations that go beyond the RPA. While none of these revelations may 
come as a surprise to the reader, it should serve as a reminder that even our most cherished approximations 
may not be controlled in any sense of the term. It is more likely that they were the first to appear in 
the literature and probably the easiest ones to use thus explaining their popularity. The dielectric function 
written down above has the attractive feature of reducing to the familiar Lindhard dielectric functionQ for 
extremely weak coupling and at the same time giving us something very different for stronger coupling. 
Furthermore, if we ask the Luttinger liquid community to write down a dielectric function of the Luttinger 
liquid, they are in all probability, going to write down the traditional RPA dielectric function'. But we 
submit that this is a serious mistake. The more nonideal a system is, the more its dielectric function differs 
from the traditional RPA, even at high density. Parenthetically, we note that having granted the fact that the 
system is nonideal (namely, no Fermi surface or close to being devoid of one) it makes no difference whether 
the system is at high or low density, the point is, it is the generalised RPA that is still valid for this system. 
The only drawback of the above approach is that if we compute the one-particle Green functions we find 
that the imaginary part vanishes identically. This is unfortunate and we have to do better in order capture 
lifetime effects. The RPA hamiltonian neglects fluctuations in the momentum distribution of the electrons. 
In order to recover a finite lifetime of single particle excitations, we find that it is important to study the 
generalised Hq rather than the more simple Hjipa- The fluctuations in the momentum distribution may be 
related to the mean by the following observation. Define, 

7V(k,k') = (nknk')-("k)(nk') (13) 
. The fluctuation in the number operator is iV(k, k) = (rt^) — (nk)^ . Since = rtk for fcrmions, we have 

7V(k,k) = (nk)(l- (nk)) (14) 

Therefore, we may conclude thaLany nonideal momentum distribution fluctuates (there are however, patho- 
logical exceptions, see footnotetj). In fact, a very nonideal momentum distribution such as one for which 
("■k) = 0.5 for most momenta has the largest fluctuation. When dealing with nonideal systems, we are 
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obliged to consider fluctuations in the momentum distribution. In order to study lifetime effects therefore, 
we have to include the full Hq. This may be written more transparently as (apart from additive constants) 
Ho = Hup A + Hfi , 



q#0 k 



k+q/2'^"-k-q/2 



(15) 



The full Fermi propagator may be evaluated by treating the fluctuation part as rS, perturbation and using the 
functional methods of Schwinger illustrated brilliantly by Kadanoff and BaymS. Before we plunge into the 
details it is important to keep in mind that the mean also changes when we consider fluctuations. That is, 



{nk)fi 



(16) 



here {n\^)RPA is given by Eq.( |^) or Eq.( pT| ) the rest is nonzero only when the fluctuation in the momentum 
distribution is large (which is, unfortunately, almost always the case when {n\^)RPA is nonideal) Therefore, 
now, (5nk refers to fluctuation around the full average rather than the RPA average. The final answers are 
given bdow, and it is hoped that the reader can rederive them using the. references quoted in the bibliography 
(mainlyOjEl). We shall adhere to the notation of Kadanoff and BaymB. In their notation the final answers 
for the single-particle Green functions are as follows (we assume in the following that F{p) ^ 0, however one 
may investigate the limit -F(p) — > 0, here z„ — {2n + 1)tt/(3), 



G„(k) - 
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and 



eic + M - S„(k) 



E„(k) = G„(k)^^(k) 



(17) 



(18) 



i^fk) 



E 

q.qVo 



y2 



7V(k q,k q) 



(19) 



From Eq.( |l8|) we may obtain the real and imaginary parts of the retarded self-energy, and from there the 
spectral function and the collision ratesj. 



Similarly, 



A{p,uj) 



Fip) 



(20) 



(21) 



if k{p,uj) ^ {lu — Ep + fi)^ — 4:F{p) < and both are zero otherwise (when F{p) ^ 0). It can be seen that 
the spectral function is peaked around ip ~ ^ with a width of the order of 2y^F{p), and the collision rate is 
vanishingly small for those values of p for which F(j3) is close to zero. It may be shown that the momentum 
distribution including possible fluctuations is given by. 
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2 f''/^ 

(np) = (-) / dd cos^e ^= 



1 



(22) 



and (lip) = for F{p) < 0. Again, it may be seen quite easily that Eq.( ^2|) is identical to Eq.( ^ when 
fluctuations in the momentum distribution are ignored. The only sticking point now is the computation of 
the fluctuation in the momentum distribrrtion. This may be done in a similar manner by employing the 
functional methods of Kadanoff and BaymB. However, we shall defer the computation of this quantity until 
the next section. Incidentally, the formulas we suggested in our preprints for this quantity in retrospect, are 
probably rather poor. Instead we shall adopt a more rigorous approach in the next section. The dielectric 
function is also modified as a result of these fluctuations and the final formula for the dielectric function 
including possible fluctuations in the momentum distribution is. 
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ee//(q,a;) = e,_«P^(q,..) - (23) 

Here, 

^ iV(k + q/2, k' + q/2) - jV(k - q/2, k' + q/2) - jV(k + q/2, k' - q/2) + iV(k - q/2, k' - q/2) 

,w ~ ek-q/2 + ek+q/2j(W - ek'-q/2 + ^k'+q/2) 



^ ^, . - -q/2 + ek+q/2 j - e^' _q/2 + Cfc' +q/2 J 



(24) 



A full derivation of this is relegated to the appendix in order to avoid interupting the flow of ideas. Assuming 
that CO has a small imaginary part, we recover both the real and imaginary parts of the full dielectric function. 

III. REPULSION-ATTRACTION DUALITY 

In this section, we point out an interesting and even a seemingly paradoxical duality between attraction 

and repulsion. In a homogeneous electron gas, electron repel each other. However, wc may show that by 
virtue of the electron being fermions, an exchange operation on the interaction part of the full hamiltonian 
leads to a change in sign of the interaction and also a change in the momentum exchanged by the interacting 
electrons. To see this, let us write down the interaction part of the full hamiltonian again. 

Hint = ^ X] 4+q/2'^k'-q/2^k'+q/2Ck-q/2 (25) 

We may rearrange the various fermion operators in this interaction so that precisely the same quantity may 
be rewritten as, 

Hint = -Yl Yl 4+q/2^I'-q/2^k'+q/2Ck-q/2 (26) 

Two important changes have occured. First there is the change in sign. This is the reason for the term 
repulsion-attraction duality. The second is the momentum that is exchanged is different. The momentum 
carried away by the virtual photon is no longer q but is k — k . The change in sign is peculiar to fermions, 
whereas the other feature is present for bosons too. Let us now use this new form and recompute the dielectric 
function. Actually, we could have done the same shuffle for the k = k part of the interaction. But we shall 
relegate a more careful examination of these issues to future publications. The main purpose of the present 
article is to lay before the reader a scheme that is robust, rich and sufficiently general so that she may apply 
the ideas to other more practical problems. In order to compute the dielectric function suggested by this 
shuffled interaction let us proceed as follows. The full hamiltonian may be written as (no(k) = nk = cj^Ck), 

H = Y^ e\no(k) - ^ ^ ^ 5no(k + q/2) 5no(k - q/2) 

k q^tO k 



- E ^E"<i(k)«-<i(i^') (27) 

k^^k' q#o 

here Hq(k) = Cj[._|_q/2''k-q/2 Let us now apply an external field, 

He,t{t) = $^(C/ext(q, t) + UU-q, i))nq(k) (28) 

k,q 

In order to proceed, we appeal to the random-phase approximation. In a suitably generalised sense, it may 
be defined to mean the following approximation : 

[nq(k),nq'(k )] = c]^+q/2Ck'-q72<^k-q/2,k'-|-q72 ~ ^k' +q72'^''~l/2'^k' -q72,k+q/2 
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« (no(k + q/2) - no(k - q/2)) 



(29) 



[nq(k),no(p)] = nq(k)(5p,k_q/2 - ^p,k+q/2) (30) 

Here we have to retain no(k) as an operator in order to include effects due to fluctuations in the momentum 
distribution. Let us now write down the equation of motion of this operator. 

«^"q(k) = (ek-q/2 - ek+q/2)"'q(k) - ^ — ^(5no(k - q/2 - q )nq(k) 

qVO 

-5no(k + q/2-q')nq(k))- ^ ^^(no(k + q/2)nq(k') - «o(k - q/2)nq(k')) 

+ (C/e.t(-q,f) + C/;,,(q,f))(no(k + q/2) - no(k- q/2)) (31) 
Write no(k) = (no(k)) + 5 no(k) 

d 

*^("q(k)) = (ek-q/2 - ek+q/2)(«q(k)) 

- E ^(('5"o(k - q/2 - q')nq(k)) - (5no(k + q/2 - q )nq(k))) 
qVo 

- E ^((^o(k + q/2)) - (no(k - q/2)))(nq(k')) 

- E ^((5no(k + q/2)nq(k'))-(5no(k-q/2)nq(k'))) 

k^ik' 

+ (t/e.t (-q, i) + C/;,, (q, f)) ((no (k + q/2)) - (no (k - q/2)) ) (32) 
d 

«^('5no(p)nq(k)) = (ek-q/2 - ek+q/2)(^no(p)"-q(k)) 
- E ^((^rio(p)^no(k - q/2 - q')) - (^no(p)feo(k + q/2 - q')))(nq(k)) 

- E ^(("o(k + q/2)) - (no(k - q/2)))(<5no(p)nq(k')) 

k#k' 

- E ^((^rio(p)5no(k + q/2))-(5no(p)5no(k-q/2)))(nq(k')) 

k^tk' 

+ (C/ext(-q,t) + C/:,t(q,t))((feo(p)5no(k + q/2)) - (,5no(p)<5no(k - q/2))) (33) 
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Assuming that, 



C^ext(-q,i) =C/ec.t(-q,0)e-'' 



(34) 



(w - ek-q/2 + ek+q/2)( ^ Mp) nq(k)) = -J2 ^(^(P,k - q/2 - q') - iV(p,k + q/2 - q'))(nq(k)) 

qVO 

-((no(k + q/2)) - (no(k- q/2))) ^ ^{S no(p)nq(k')) 
-(7V(p,k+q/2)-iV(p,k-q/2)) ^ ^-^{n^(k')) 

kVk 

+ [/ext (-q, 0)(Ar(p, k + q/2) - N{p, k - q/2)) (35) 
(w - ek-q/2 + ek+q/2)(nq(k)) = - ^ ^((^no(k - q/2 - q')nq(k)) - (5no(k + q/2 - q')nq(k))) 

qVO 

-((no(k + q/2))-(no(k-q/2))) ^ !^(n<,(k')) 

kVk 

- E ^((^«o(k + q/2)nq(k')) - {5no{k - q/2)n^{k ))) 

kVk 

+ f/ext(-q,0)((no(k + q/2)) - (no(k-q/2))) (36) 

In order to simplify this further, let us assume that we are in the weakly nonideal regime. That is, it is 
legitimate to treat the momentum distribution as possessing a sharp Fermi surface and no other striking 
features. Let us now define. 



S no(p) nq(r)) = -1^ ^-f f ^'r' v{r') e*-(^^'-^^(7V(p,k-q/2-q')-iV(p,k+q/2-q')) {h^{7^)) 

/ E^"''^'''"'^(<"o(k + q/2))-(no(k-q/2)))t;(r')(5no(p)nq(r-')) 
k 

k 

+ f/ext (-q, 0) i ^ e-'''-"^(iV(p, k + q/2) - 7V(p, k - q/2)) (38) 

k 
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+ i^){fi^{f)) ^e*-(^^'-^^)((feo(k - q/2 - q>q(k)) - (<5no(k + q/2 - q )nq(f'))) 

k,qVO 

-^E / ^(r-')e*-(^^'-^^((no(k + q/2))-(no(k-q/2)))K(f')) 

_l_J2^u.ir'-f) f ^3/ ^(r')((5no(k + q/2)nq(r'))-(Jno(k-q/2)nq(r'))) 
k 

+ f/,,,(-q,0)-i ^ e-''^-^^((no(k + q/2)) - (no(k - q/2))) (39) 



V 

k 



Let us now compute the following quantity, 

A(^) = ^ E e'^-^'-'^iinoO^ + q/2)) - (no(k - q/2))) (40) 

k 

where R = r — r. If A:/ is sufficiently large ( kf » q, kf » 1/R) we may write, 

U{R) « Coci.Vr6^R) (41) 

Co = ^^C-)l^dR ^(sinikfR) - {kfR)cos{ksR)){sm{^) - {^)cos{^)) (42) 
(a;+i^)(5no(p) h^{r)) = ~ I ^(^') e^'^-^'''-'^(^(p, k-q/2-q')-A^(P, k+q/2-q')) (nq(r')) 

+Coq.Vf^(r)(5no(p)nq(r)) 

Z'^'^' E«"''^'''"'^(^(P''^ + q/2)-^(P'k-q/2)Mr-')(nq(r)) 
k 

+ t/ext(-q, 0) :^ ^ e-*-^"(iV(p, k + q/2) - iV(p, k - q/2)) (43) 

k 

(a;+i^)(fiq(r-))=-i ^ /c^V' t;(r>*-(^^'-^^((^no(k-q/2-q')nq(r'))-(^no(k+q/2-q')nq(r'))) 

k.qVO 

+Coq.Vf^(r)(fiq(r)) 

_ 1 ^g*.(f'-fO /■ ^3/ „(/)((5no(k + q/2)nq(r-')) -(5no(k-q/2)nq(r-'))) 
k 

+ C/e.t(-q, 0)i ^ e-^'^-^"((no(k + q/2)> - (no(k - q/2))) (44) 

k 

At this stage the author is unable to complete this calculation, but hopefully the reader appreciates the 
spirit of the discussion. The author apologises in advance for this omission. 
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IV. IMPROVED THEORY OF SINGLE-PARTICLE PROPERTIES USING BOSONIZATION 



In this section, we use the sea-boson method suitably generahsed to accomodate fluctuations in the momen- 
tum distribution to compute single-particle properties. In the appendix we show how to derive the diclcctpc 
function using this method. The dielectric function may be evaluated by more conventional means as wellllil. 
However, the single-particle properties evaluated using the sea-boson method is more accurate as one is 
able to treat the problem more systematically and the physical interpretation of the formulas is also more 
transparent. The really systematic approach would be to write down the sea-boson correspondence suitably 
generalised to accomodate fluctuations in the momentum distribution and evaluate the various propagators 
like we did in our earlier worka. However, we shall adopt a simpler but hopefully the correct approach here. 
The idea is to simply borrow from our earlier work except that the zeros of the dielectric function are no 
longer the zeros of the RPA dielectric function but that of eg/ / and we need to interpret the derivate of the 
dielectric funcion with respect to frequency that appears in these formulas as being the derivative of the full 
dielectric function. Therefore, let us write down the various formulas. The full hamiltonian in diagonalised 
form has the following appearence, 

=^c^,(q)4(q)d,(q) (45) 
nk = n'^(k) i^i(k) + (1 - n^(k)) i^2(k) (46) 



where. 



^iW = : ^ (47) 



1 



l+S^(k) 



^^C') - V-TTs^ (48) 



^ (wi(-q) +k.q/m-eq)2 
(t^i(-q) + k.q/m + eq)^ 



k 



It may be observed that this quantity is just the frequency derivative of the polarization. In other words, 

V d 

9i'^{^) = — (^).^=a;,(q)<-flPA(q>^) (52) 

By analogy we may generalise this so that when fluctuations are introduced, 

V d 

ff^"'(q) = —{-Q^)o.=u.M)^efM,^) (53) 

We shall employ this latter definition in our analysis. Wi(q) is the zero of the overall dielectric function, 

ef^^(q,c.,)-0 (54) 
It is worthwhile simplifying the effective dielectric function. 



9 



eeff{q,uj) ^ eg-RPA{q,i^) + — '—7- — ^ (55) 



The principal part of this is given by, 



«-p,PA(q>^) - i)«-i?PA(q>^) - ^^{(i,^)) 

e;//iq,t^j = Cg-flPAiq^^J + 

where, 



ee//(q, ^) = fg-B.PAi(l, ^) + jp 7—^ (56) 



Eg-PPA (q, t^) = 1 + 2^ '- ^ '- (57) 



_ /3 

ep{q,Lu) = l + f}^ (58) 

k ra 

and the principal parts of these functions are just the real parts. 

er-pP.(q, ^) = 1 + ^ E ^ %+q/.J^.-q/. (59) 

k ^ m 



A. Just How Many Zeros Does the Dielectric Function Really Have ? 

In this section, we address what is perhaps the most vexing problem in this approach. What is the size of 
the set to which i in Eq. ( ^) belongs ? If one counts only the collective mode which is what a naive approach 
would lead us to do, then we are ignoring a large (infinite !) number of particle-hole modes. And yet, it 
seems that the particle-hole modes don't come about naturally and must be forced into the formalismu. It 
is really important to address this issue since ignoring it would mean that the excitation spectrum of the 
homogeneous electron gas possesses a gap when in fact it should not. In other words, if one counts only 
the collective mode(plasmon) one arrives at the unavoidable conclusion that the excitation spectrum has a 
finite gap(equal to the plasmon energy). How then should one introduce the particle-hole mode ? In this 
section, we show how to introduce both the particle-hole mode and collective mode in a unified manner. In 
our previous worka we suggested that the zeros of the RPA dielectric function should be interpreted as the 
maxima of the dynamical structure factor. This is quite an unusual and drastic departure from the notion 
of a root of a function. Being physicists we accept it as it has a physical interpretation. Let us examine the 
diagonalised form of the full hamiltonianQ: 

-^^,(q)4(q)d,(q) (60) 

For small q, we expect this sum to involve just the collective mode. In other words, 

Hfuu^ E '^c(q)4(q)4(q) (61) 

c{— small 

This may be achieved by employing the following device, 

Hfuu = E / "^^^ W{(i,iu) 4(q)^^-(q) (62) 

Now these operators, the dressed sea-bosons in the new language obey a different kind of commutation rule: 

[d^(q),d^,(q')]=0 (63) 
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K(q),dl,(q)]=<5^^„'<5q,q' (64) 

where 6^ j is a Kronecker delta function. Also the object VF(q, oj) is a weight suitably chosen so that the 
the sum reduces to just the one over the collective mode in the small q limit. It is also clear that M^(q, w) 
has to have dimensions of inverse energy. Let us therefore postulate a form, 

W^(q,c^) = -co/m(^^) (65) 
where cq is a suitable constant. We know from textbooks! that(in 3D), 

Limitq-^QW{(\Tijj) — —Co Limitq-^o r) = ° ^ 5[ijJ — ujp) (66) 

e(q, w) 2 

In order that we reproduce the right collective mode we must choose, 

CO = — (67) 



P 

The choices for W^(q, lu) in other dimensions are worked out in the appendix. Now we have in our hands a 
convenient way of labeling excited states of our system. The eigenstates of the system have energy eigenvalues 
given by, 

n^{q) ^LuW{q,uj) Alo (68) 

The spacing Alo is the smallest possible spacing between the energies, which is arbitrarily small. Thus we 
can see that, in most cases it does not cost a finite energy to excite the system. This corresponds to the 
particle-hole mode. On the other hand, if the weight W^(q, lu) is singular as it happens in the vicinity of 
oj = (jjpand q = 0, we have a finite gap. 

^uj^ujp (q w 0) (5(w - ujp) Alu (69) 

Now if Lim^^ajp^i'-^ ~ ^p) — 1 (this defines Auj if you like), then we see the emergence of the collective 
mode. In gapped systems we expect only the collective mode and no particle hole mode. This is possible 
only if for any lj ^ Q, 

H^(q,cj)=oo (70) 

For the homogeneous Fermi system, this situation corresponds to Wigner crystallisation. In a Wigner crystal, 
in order to create excited states, you need to create phonons which require a nonzero-amount of energy each 
time. 

Thus, all energies are allowed but each comes with a "weight" corresponding to how strong the structure 
factor is at that energy. It may be seen that in 3-dimensions this would mean that for small q the sum 
over i is just the collective mode, but for larger q, we start summing the particle hole modes as well. The 
sums in the objects Sa and Sb have to be reinterpreted as well. It is rather unfortunate that simple and 
straightforward interpretations are not possible, and one is forced to take recourse to devious means such as 
the one we shall now describe. Let us postulate, 

/"OO - 

5.(k).^| ^Wi^..)j-^±l^,fM) (71) 

5b (k) =Y, dcj W{q,c.) ., .9^(q) (72) 

a'M) = "^^^i — 7 (73) 



11 



and, 

poo 

W^(q,w) =VK(q,cj)/ / dui W{q,Lu) (74) 
Jo 

This seemingly adhoc ansatz becomes more plausible when one realises that in the small q limit it has a 
familiar form corresponding to the collective mode. 

ci— small ^ 
q— small ^ 

Let US now move on to the full propagator. 

B. The Full Propagator 

Let US borrow from our earlier workB and write down the final formulas, 

(77) 
(78) 

In the above formula, the index i runs over both the collective mode as well as the particle-hole modes(i = 
c, ki). As we have pointed out earlier, this must be reinterpreted to mean a weighted sum with the dynamical 
structure factor suitably normalised appearing as the weight. 

<?k.q(x) = ~e-^ <i-x(_i_)AO (-q)c.k(q) + ^ C/q(x)AO (-q) = -/^,q(x) (79) 



Set i7q(x) =e-* '^■"t/o(q) 



».M = (ji7T)E(^)'(A£(-q))' (81) 

1 k 



as 



^2(q)-(^)E(^k(-q))' (82) 

k 

Here A^(— q) is the A with a nonintcracting momentum distribution. The interacting coefficients are given 
t/k,q(x) = /^,q(x)[ak(q),6l(q)] + /k,-q(x)[ak(-q), 6.(q)] (83) 
7^o = ea:p(- ^ /k^q(x)/k,q(x)[6j(q), a]^(q)][ak(q), 6|(q)]) 

k,q,i 
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xexpi-l E /k,q(x)/k,-q(x)[ak(-q),6i(q)][ak(q),6l(q)]) 



2 

k,q,i 



X ea;p(-J ^ /k,q(x)/k,-q(x)[ak(-q), &i(q)][ak(q), 6l(q)]) (84) 



2 

k,q,i 



The commutators are given as before, 



[6,(q),aUq)] = (-M4^)g,(q) = [ak(q), 6l(q)] 



(86) 



Mq) 

^Wi(q) 



[6i(q),ak(-q)] = -( , . k.q )gi(q) (87) 



[ak(q),6,(-q)] = (-AL^)5,(-q) (88) 



5?(q) = ^^T:r-T (89) 



£|c.=a;,er//(q,w) 



C. Exchange Effects within RPA 

It may be possible to employ a myriad of other similar approaches each differing from the other in how 
the notion of RPA is implemented. To give an example, let us rewrite the full hamiltonian differently. 



H = 2^ CkCfcCk - 2^ 2^ nk+q/2nk-q/2 



k q^tO 

+ E^EMk)n-q(k') (90) 

q#0 k^tk' 

Here the exchange term has been singled out for special treatment. We could rewrite the last term differently, 

rr Y^ t Y^ ^(q) Y^ 

H = }_^ CkCj^Ck - 2^ 2^ nk+q/2nk-q/2 

k qytO k 

- E ^ E "k-k'(k/2 + k72 + q/2)nk'_k(k/2 + ],/2- q/2) (91) 

qT^O k^tk' 

This may be interpreted as an exchange term. The negative sign suggests precisely this. These two hamil- 
tonians lead to different answers when the RPA is carried out on them. In fact, even the exchange term 
involving just the number operators can be treated in two different ways. One way is to retain it as it is, 
the other is to rewrite it as. 
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Hex = ^k4ck - Y ^ S nk-q/2 (92) 

k q^^O k 

The effective dispersion includes the exchange energy, 

?k = ek-5]^K-q) (93) 
q#0 

and S rik — ni^ — fik Again when the RPA is carried out on them they yield different answers. One has to 
examine this issue more closely. In particular the following questions spring to mind. Should we use one or 
the other or a combination of the two ? The answer may he .given by appealing to our bosonization scheme. 
For more details, the reader is refered to our earlier worksH ,t3. Just to make this article self-contained, it is 
desirable to fix some notation. 

no(k) = c[ck (94) 

the ak(q) is a sea-boson corresponding to the fermions Ck- The coefficients Ak(q) — \/'T-k+q/2(l ^ ^'ik-q/2) 
are c-numbers for now. 

4+q/2Ck-q/2 = Ak(q)ak(-q) + aj^(q)Ak(-q) (95) 

Let us now try and figure out how to decompose products of four fermions, (let us further assume that none 
of the four indices are equal to any other index) (q 7^ and k 7^ k ) 

4+q/2'=k-q/2cj^'_q/2'^k'+q/2 = I^k (q)ak (-q) + aj, (q) Ak (-q)] [A^' (-q)ak' (q) + a|^, (-q) A^' (q)] (96) 

the above equality is within RPA. It may be noted that a slightly different regrouping yields a totally different 
formula, 

4+q/2'^k-q/24'_q/2"^k'+q/2 = ^4+q/2^k' +q/2Ck' _q/2'^k-q/2 

= -[Ak k' , (k-k')a^ „(k'-k) + a'f , (k-k')A^ Jk' - k)] 
^ i+>|-+q/2^ ^ + ' i+^+i^ 

x[A^^_„(k'-k)a^^_„(k-k') + at , Jk'-k)A^^ Jk-k')] (97) 

2^22 2'22 ■2+"2 T 2'22 

Let us now examine the commutator obtained using Fermi algebra, 

[4-q/2Ck+q/2' 4+q/2^k-q/2<'_q/2Ck'+q/2] = K (k - q/2) - riQ (k + q/2)] cj^, C^' +q/2 

+ [<5k'-q/2,k+q/2(l " "o(k - q/2))4+q/2Ck'+q/2 + '5k-q/2,k'+q/2"o(k -|- q/2)c|^, _^^2Ck-q/2] (98) 

It is clear that if one uses the ansatz in Eq. ( |9^) then one obtains only the first term in Eq. ( |9|) . In order to 
obtain the second term it is necessary to go beyond this and try and include some portion of Eq. ( |9^) . Let 
us first try to use the ansatz in Eq. ( ^) while completely ignoring the ansatz of Eq. ( ^) In this case the 
interaction term takes a different form, 

ff/ = -E E ^[Ak(q)ak(-q) + aUq)Ak(-q)][Ak'(-q)ak'(q)+at,(-q)Ak,(q)] (99) 
q#o k^iik' 

Here two important changes have occured. First, the interaction Wq is replaced by Wk_k' there is a 
change in sign. The physical meaning of these mysterious changes will be defered until later when we have 
investigated the details. Let us now write down the full hamiltonian. 
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H = J2^^i<l)4i<l)M^) -EE ^[Ak(q)ak(-q) +at(q)Ak(-q)][Ak'(-q)ak'(q) +at,(-q)Ak.(q)] 

(100) 

Here, 

'ik(q) = ek+q/2 - ek-q/2 (101) 

= ek - E T7^k-q (102) 

q#0 

Let us diagonalise this via a Bogoliubov transformation. To this end one assumes that the diagonalised form 
has the following apperance, 

F = ^a;i(q)4(q)rfi(q) (103) 

i,q 

Wi(q)(ii(q) = E'^k(q)[di(q),ak(q)]ak(q) +E'^k(-^)['^i(^)'"k(-q)]ak(-<l) 

k,q k,q 

-EE ^(Ak'(q)Mi(q),«k'(-q)] + ['^i(q),<'(q)]Ak'(-q))(Ak(-q)ak(q) + 4(-q)Ak(q)) (104) 
q^^o k^^k' 

'^i(q)[c;i(q),ak(q)] = '^k(q)[rfi(q),oj^(q)] 

- Ak(-q) J2 (q)Mi(q), (-q)] + [^.(q), 4' (q)]Ak' (-q)) (io5) 

'^i(q)Mi(q),ak(-q)] = -a)k(-q)[di(q),ak(-q)] 
+ Ak(q) E ^(Ak'(q)[rfi(q),«k'(-q)] + [rfi(q),at,(q)]Ak'(-q)) (106) 

[dM), 4 (q)] = .^^^~^\ , [ df e''^-^^;(r){i?i(q, f) + R^id, r)} (107) 

Mi(q), ak(-q)] = , , [ df e*-~z;(r){i?i(q, r) + i?2(q, r)} (108) 

Wi(qj +a;k(-qj J 

Ri{<i,r) = :^E^""''-''V(q)['^i(q),«k'(-q)] (io9) 

k' 

^2(q,r-) = :^^e-''''-X'(-q)[rfi(q),4'(q)] (^0) 



Here, 



V , 

k 
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Now define, 



(a;,(q) -(i,v^_, (q))i?2(q,r ) = -J df i^i(f- f' ; q>(f){i?i(q, f) + i?2(q, r)} (111) 
(a;i(q)-a;iv^_,(q))i?i(qy)= J rfr F2(r - r'; q)t;(f){i?i(q, f) + i?2(q, r)} (112) 
Fi(f - r'sq) = 1 ^ Ak(-q)e*-(^^-^^') (113) 

k 

F2(r - r'sq) = 1 ^ Ak(q)e*-(" > (114) 

k 

i?i(q,f) + i?2(q,r') = i?(q,r) (115) 
{tOiid) - uji^^ (q))i?(q,f' ) = j d\ /o(q;r'- f>(r*)i?(q, f) (116) 

/o(q;r - r') = 1 ^ n^(k)e'"^-(^^-^')(-2i)sm(iq.(r - r')) (117) 

{^ny \r — r\ \r — r\ \r — r 

(^.(q) _ iS:Zfl)ij(q,f') = _i Coici){ci.V^)v{r^)Rici,f^) (119) 

[a;,(q) - i^:^ + i C7o(q)(q.V,^^;(r'))]i?(q, r') = -i Co(q)t;(r'')(q.V,- )i?(q, r') (120) 

This eigenvalue problem admits many eigenvalues. Actually, if not for the constraint of the square integra- 
bility of -R(q, r), for each q any positive real number would be an eigenvalue. 

[Ui{ci) - iS:^ +i Co(q)(q.V-t;(f'))]i?(q,f') = -i Co(q)t;(f')(q.V-)E(q, f') (121) 
Let us now solve this equation, 

R{ci,f)=R{ci,fo)exp{-. / dr , ) (122) 

iqcosej^^ l--Coiq)v{r) 

For some suitable ro- 

This work was supported by the Dept. of Physics at University of Illinois at Urbana- Champaign. The 
authors may be contacted at the e-mail address setlur@mrlxpa.mrl.uiuc.edu. 
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V. APPENDIX 



In this appendix, we show how to derive the full dielectric function using the generalised RPA. Along the 
way we point out some pitfalls and possible generalisations. Let us write the generalised RPA hamiltonian 
and try and compute the dielectric function. 

Ho = Yl ^kno(k) -J2^J2M^ + q/2)no(k - q/2) (123) 

k qT^O k 

He.t{t) = E(t^ext (qt) + UU-qt)) ^ nq(k) (124) 
q k 

here nq(k) = c^_^_q/2^i'fk-ci/2 and no(k) = c|^Ck Let us now write down the quation of motion for nq(k). 

d I 
^^"qW = X!^k'K(k),no(k )] 

k' 

-E^EMk),no(k')]no(k'-q') 

qVO k' 

+ E (Ue.t{cit) + C/:,,(-q t))[nq(k), n^, (k')] (125) 

k'.q' 



Now, 



K(k),nq'(k')] = [4+q/2Ck-q/2,4'+ '/2Ck'-q'/2] 



- 4+q/2Ck'-q72^k-q/2,k'+q72 ~ c|^'+q'/2'^k-q/2'5k+q/2,k'-q72 (126) 

One may approximate this as, 

[nq(k),nq-(k')] = [4+q/2Ck-q/2,cj^'+q'/2Ck'-q72] 

= [no(k + q/2) - no(k - q/2)]5k,k'^q,-q' (127) 

The next approximation would be to replace the number operator by its c-number expectation value. But 
we shall desist from that for the moment. 

K(k),no(k')] = nq(k)(5k',k-q/2 - 4',k+q/2) (128) 



(k) = (ek-q/2 - ek+q/2)"-q(k) - ^ -^(jiQ^k. - q/2 - q ) - no(k + q/2 - q ))nq(k) 

qVo 



+ (I7e.t(-q0 + J7*,,(qt))(no(k + q/2) - no(k - q/2)) (129) 

(nq(k)) = C/e,t(-qt)Cq(k) + U:^,{c^)D^{\^) (130) 
Let us now ignore fluctuations in the momentum distribution. That is, we are allowed to replace 
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(no(k')nq(k)) = (no(k'))(nq(k)) 



(131) 



w f/ext(-qt)Cq(k)-a; [/;^t(qi)£)q(k) = (ek-q/2-ek+q/2)t^ext(-qt)C'q(k) + (ek-q/2-ek+q/2)f/ext(q*)^q(k) 
- E ^(("o(k - q/2 - q')) - (no(k + q/2 - q')))[/e,t(-qt)Cq(k) 

- E ^(("o(k - q/2 - q )) - (no(k + q/2 - q )))C/;,,(qt)Dq(k) 

qVO 

+ (f/ext(-qi) + [/e*.t(qt))((rio(k + q/2)) - (no(k - q/2))) (132) 
w Cq(k) = (ek-q/2 - ek+q/2)Cq(k) - ^ ^((no(k - q/2 - q')) - (no(k + q/2 - q')))Cq(k) 



Since, 



+ ((no(k + q/2)) - (no(k - q/2))) (133) 
UefMt) = Ue,t{cit) + ^(p'_q)f/ext(qt) (134) 

(/''-q) = E^-^('^) (135) 



u - ek+q/2 + ek-q/2 

If we use the definition of the dielectric function we get a wrong answer. 



(136) 



, X C^ext(qi) 1 /1q-7^ 

eWRONGici, OJ) = = (no(k-q/2))-(no(k+q/2)) 

eJJVil i-l-yZ^k w-ek+„/2+?k-<,/2 

The reason is because we have been very cavalier in our treatment of correlations. One has to consider the 
full Hamiltonian when dealing with the dielectric function rather than just Hq. Let us now try and do this. 

»^<(k) = (ek-q/2 - ek+q/2)<(k) + ^ ^[«q(k), Pq']p*_q- + ^<'K(k)>P-q'] 

qVO qVO 

+ E (^e.t(q'i) + C/:.t(-q't))K(k), Pq'] (138) 

qVO 



["q(k),Pq'] = E['=k+q/2Ck-q/2,4'+q'/2Ck'_q'/2] 



k' 



E^k+q/2^k'-q72'^k-q/2,k'+q72 ~ E 4'+q72'^'*-'l/2^k+q/2,k'-q72 
k' k' 
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~ '5q,-q' ["o(k + q/2) - no(k - q/2)] 



(139) 



= (ek-q/2 - ek+q/2)<(k) + y (no(k + q/2) - no(k - q/2))p^ 

+ («7ext(-qi) + [/e*c.t(qt))(«o(k + q/2) - no(k - q/2)) (140) 

Let us make a first pass at the computation of the dielectric function. Here, we make use of mean-field 
theory, that is, replace 

(no(k')pq) = (no(k'))(pq) (141) 
u; «(k)) = (ek_q/2 - ek+q/2)«(k)) + ^((no(k + q/2)) - (no(k - q/2)))(p^) 

+ (f/ext(-qi) + [/ext(qt))((rio(k + q/2)) - (no(k - q/2))) (142) 



t ^^ _ ^q (no(k + q/2)) - (no(k - q/2)) ^ 

V Ul- ek-q/2 + ek+q/2 



(t/..,(-.<) + t/;..M) '"°"' + -'/^)'-<"°"'-/^'> (143, 

- ^k-q/2 + ek+q/2 



This means, 



Po(q,^) 

e(q,w) 



(p-q) = C/e.t(qi)^}e^ (144) 



Po(q,a;)=V^"°('-q/^)^-^--°('- + q/^)^ (145) 

■V^ W - ek+q/2 + ek-q/2 



e(q,w) = l-^Po(q,w) (146) 



From this and the fact that 



/ N C^ext(qi) , n 

eg-flPA(q, w) = ; = e(q, w) (147) 
t/e//(qt) 

Next we would like to include fluctuations. Let us do this differently this time via the use of the BBGKY 
heirarchy. 

i-«(k)) = (ek-q/2 - ek+q/2)«(k)) + ^((no(k + q/2)) - (no(k - q/2)))(p^) 

+ ^(P2A(k + q/2, q) - F2A(k - q/2, q)) 
+ (C/ext(-qi) + [/ex((qt))(no(k + q/2) - no(k - q/2)) (148) 

Here, 

F2A(k',q;t) = (no(k')p^) - (no(k'))(pJj) (149) 
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F2(k';k,q;i) = (no(k')n^(k)) - (no(k'))(n^(k)) (150) 
J^2A(k',q;t)=^F2(k';k,q;t) (151) 



i-F2(k';k,q;t) = (ek_q/2 - ek+q/2)^2(k';k,q;t) + ^ (7V(k' , k + q/2) - iV(k',k - q/2))(p^) 



V, 



+ -^((no(k + q/2)) - (no(k- q/2)))F2A(k ,q;t) + (f/ext(-qi) + U:,,{cit)){N{k ,k + q/2) - Ar(k ,k- q/2)) 

(152) 



Let us write, 



J^2(k';k,q;t) = t/e,t(-qt)i^2,a(k';k,q) + C/;,,(qt)J^2,6(k';k,q) (153) 



Also define, 



7V(k,k ) = (no(k)no(k )) - (no(k))(no(k )) (155) 
CO F2,„(k';k,q) = (ek-q/2 - ek+q/2)F2,a(k'; k, q) + ^(Ar(k',k + q/2) - Ar(k',k - q/2))(pq) 

+ ^((no(k + q/2)) - (no(k - q/2)))F2A(k', q) + (A^(k', k + q/2) - iV(k', k - q/2)) (156) 

-oj J^2,6(k';k,q) = (ek_q/2 - ek+q/2)i^2,6(k'; k, q) + ^ (7V(k' , k + q/2) - Ar(k', k - q/2))(p;) 



+ ^((no(k + q/2)) - (no(k - q/2)))F|^(k , q) + (iV(k , k + q/2) - N(k , k - q/2)) (157) 

f(a u.) F» (k V ^(k',k + q/2)-jV(k',k-q/2) , 
e(q,w) i'2A(k ,q) - t7 2^ ~ — Z ^P^' 

^ ^ ek-q/2 + ek+q/2 

^ jV(k\k + q/2)-A-(k\k-q/2) ^^gg^ 

^ W - ek_q/2 + ek+q/2 

o; Cq(k) = (ek_q/2-ek+q/2)Cq(k) + ^((no(k+q/2))-(no(k-q/2)))(p;) + ^(F£4(k+q/2)-J^2"A(k-q/2)) 

+ ((no(k + q/2))-(no(k-q/2))) (159) 
r _ ^;q (no(k + q/2))-(no(k-q/2)) ^^. 

V U - ek_q/2 + ek+q/2 
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(no(k + q/2))-(no(k~q/2)) 

^ ^ £k-q/2 + ^k+q/2 



^ J-°^(k + q/2, q) - f °^(k ^ q/2, q) ^^^^^ 



— ek-q/2 + ek+q/2 



After all this, it may be shown that the overall dielectric function including possible fluctutations in the 
momentum distribution is given by, 

e,,,(q,c.) = e,_HP^(q,.) - (^)^ ^'^""f . (161) 

Here, 

A^(k + q/2, k' + q/2) - A^(k - q/2, k' + q/2) - A^(k + q/2, k' - q/2) + A^(k - q/2, k' - q/2) 

- '^k-q/2 + ek+q/2)(t^ - V -q/2 + V+q/2) 



P2(q,c.) = ^ 



(162) 



eg-flPA(q, ^) - 1 + TF 2^ — (163) 



k 



£k+q/2 + £k-q/2 
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The pathological exception is for nonideal distributions such as, (rik) = for |k| < 0.5 fc/ and (rik) = 1 for 
0.5 fc/ < [kj < 1.5 kf and (rik) = for |k| > 1.5 fc/ this is very nonideal but has zero fluctuation, but hopefully 
these are not physical. 
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